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The Raman scattering (RS) spectra of rhombohedral LaMnOs, is considered. It is shown that 
broad high intensity RS line at 640 ana"^ IJ could be explained in the framework of cooperative 
, dynamical Jahn- Teller effect. This hypothesis is illustrated through the lattice dynamics calculations 

. explicitly allowing for the electron-phonon term both in crystal energy and in dynamical matrix. 

' The averaging scheme for the one mode cooperative dynamical effect is proposed 
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INTRODUCTION 



The colossal magnetoresistance (CMR) is the main feature of lanthanum manganites 0. Because of metal to insu- 
lator transition, the doped rhombohedral manganites are of grate interest. The nature of ground state in manganites 
is still under the question. There are strong coupling between orbital (electronic) and lattice degrees of freedom in 
. manganites due to Jahn- Teller (JT) effect on Mn^+ ions. Such coupling gives rise to the formation of vibronic ground 
state Isj . The vibronic ground state fundamentally transforms magnetic and transport properties of manganites 



^ ' The pure LaMnOz is the model compound for investigating cooperative JT effect. Undoped LaMnOz crystallizes 
Q in two phases: orthorhombic and rhombohedral The orthorhombic one features by the static cooperative JT 

^ ' effect. The static cooperative JT effect and i?3c space symmetry of rhombohedral phase are incompatible because 
the JT distortions break this symmetry. The only way to conserve the i?3c space symmetry in rhombohedral phase 
is dynamical JT effect. The possibility of dynamical JT effect in manganites is widely discussed 0, 0, 0- The 
^ . dynamical JT effect manifests in two basic forms: uncorrclatcd local dynamical JT distortions of [MnOg] complexes 
' and cooperative dynamical JT distortions, which can be called cooperative dynamical JT effect The most 

. powerful tool that could determine presence of low symmetry structure is Raman scattering (RS). We suppose that 
rhombohedral LaMnOj exhibits cooperative dynamical JT effect. We show that anomalous intensity of RS lines 
: observed in rhombohedral LaMnOs pj,lll| could be explained in terms of cooperative dynamical JT effect. For the 
^f) ' lattice dynamics calculation the interionic pair potential shell model approximation with electron-phonon (JT) term, 
, which explicitly allows for the dynamical effect, is used [l^ H^. 



SYMMETRY AND QUALITATIVE ANALYSIS OF RHOMBOHEDRAL LaMnOs 



B] 

' ^ As was mentioned by Iliev and coworkers P, [3| the Raman spectra of rhombohedral LaMnO^ differs strongly from 

O RS of non JT LaAlOs 

O From the symmetry analysis we know that in the i?3c space group only five modes are active in RS {Aig + ^Eg). 

^> . The experiment on LaAlO^, Q clearly shows 4 Raman active mode and one, that has very low frequency of 33 cm^^ 
was also observed Comparing the RS of rhombohedral LaMnO^ one should note strong distinction to RS of 

LaAlO^. The main difference is appearance of new intensive high-frequency line near 640 cm^^ and relative increase 
' of intensities of high frequency lines. The appearance of new RS line is clear evidence for lowering symmetry, i.e. 
that the real space symmetry of rhombohedral LaMnO^ is lower than i?3c. The symmetry lowering could be due 
to cooperative JT effect. The only JT active mode in F point of Brillouin zone (BZ) of i?3c space group is double 
degenerate Eg mode, which corresponds to the {fciajrs mode of ideal perovskite (in Kovalev notation [l^ ') or R12 (in 
BSW notation). Two components of Eg mode could be indexed as qg and (je (Fig. The symmetry group of new 
phase depends on the ratio of qg to q^ . This ratio is convenient to describe by an $ angle, defined as qg = p cos $ and 
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qe = psin$. If $ = ^, ^ ± ^ the structure will be with P2i/h space symmetry. In other case space symmetry will 
be PT. 

The JT distortions in P2i/h space group differs strongly from that in orthorhombic phase with Pnma space group: 
in orthorhombic phase the Qe distortions correspond to the M point of ideal perovskite BZ ({fciijrs in Kovalev 
notation [l^l. On other hand, in experimentally observed monoclinic phase structure 0/ JT distortions correspond 
to the Ri2 mode of ideal perovskite. 

P2i/b and PI groups both have all 12 gerade modes active in RS. These modes correspond to the following i?3c 
vibrations: split four Eg double degenerate mode; one Aig mode, and three modes. The symmetry breaking 
would give rise to the appearance of a new line. It corresponds to the A2g vibration. The intensities calculation 
revealed the remaining two A2g modes to be low RS active. For the investigation of such phase transition we use 
lattice structure and dynamics calculations. 



LATTICE STRUCTURE AND DYNAMICS CALCULATIONS 



Energy calculation model 

In lattice structure calculation the crystal lattice parameters are obtained by the optimization of crystal energy 
with respect to displacements of ions that do not change crystal lattice symmetry. The crystal energy is calculated in 
the framework of ionic pair potential shell model explicitly allowing for the JT term, which is simulated as the sum 
of lower branches of adiabatic potentials of single [MnOe] complexes 0, 

E = ELat + EjT (1) 

In the pairwise-potential approximation and the shell model, the lattice energy can be written as 

where Ecoui describes Coulomb interaction in crystal and the index i specifies ions in the primitive cell, the index 
k enumerates all ions of the crystal. The term kiSf describes the interaction energy between the core of the ith ion 
and its shell shifted relative to the core by an amount of 5i. Short range potential is represented by the Born-Mayer, 
screening and Van-der-Waals terms 

Vik = Cik e^]i{-Dik\rik - 5j + 4|) (3) 

Aik , „ s Aifc 

exp(-B^fcrjfc) g- 

rik f^fc 

The screening term describes the screening of Coulomb potential due to overlap of electronic densities of ion. Pa- 
rameters of pair potential model are listed in table ^ The many-particle JT contribution to the crystal energy is 
approximated by the sum of expressions for the lower branches of adiabatic potentials of the [Af nOg] clusters 



EjT = -\Ve\Y.^Ql{n)^Q'i(n) (4) 

n 

where is the linear JT interaction constant (|Ve| — 1.3 eV/A 0|), and Qein), Qe{n) are projection of displacements 
of oxygen ions onto local JT distortions. Actually, these local distortions are projection of positions of all ions 
onto [MnOe] cluster local normal JT active modes. The operation of projection can be performed using Van-Fleck 
coefficients 

where n - indexes JT clusters, C^^ e)„fc^ " is the projection of displacements caused by crystal vibration with wave 
vector k and number 7 onto local Q(e,e) distortion. 
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Raman scattering spectrum calculation model 

To calculate the phonon frequencies, one should diagonalize the dynamic matrix of the crystal. 

D = M (f^^ - F^'^ [F^'^y^ F'^^ + F-^'^) (6) 



where Mka,ni3 — Sa/slmkmn]^^/^; F'-^'^ and F^^ are the matrices of the second derivatives of the crystal energy 
(without the JT contribution) with respect to the displacements of the cores and shells of ions, respectively; F'^^ and 
F"^*^ are the matrices of the mixed derivatives of the crystal energy (without the JT contribution) with respect to the 
displacements of the cores and shells of ions; and F"'^ is the matrix of the second derivatives of the JT contribution 
in Eq. ^with respect to the ion core displacements. 

The ions shells simulate the polarizability of ions, which arises due to transitions to the excited states with different 
parity. The RS process involves virtual transitions to these states. That is why the RS intensities could be estimated 
in shell model. Hence, the intensity can be approximated concerning the crystal like classical media. 

/.ocf^V (7) 



dq 

where Iq is the intensity of mode q, aoc is the high frequency polarizability. The derivative is taken with respect to 
the amplitude of vibration q. In the shell model approximation the derivative could be calculated as follows: 




where f2 is the volume of primitive cell, U is the matrix which diagonalizes dynamical matrix of the crystal (Eq. ^ , i 
- indexes vibrational modes, a and P index Cartesian coordinates, k and m index ions in the primitive cell, Yk - is the 
charge of the shell with number fc, nik is the mass of the ion with number k. F^'-'^^ is the matrix of third derivatives 
taken with respect to the shell {S) and core (C) displacements. 

Using the same approximation for calculation of frequencies and RS lines intensities we can perform the assignment 
of RS peaks to the calculated phonons. Moreover, this approximation allows investigating the RS lines intensities 
dependence on the values model parameters like pair potentials or JT term parameters as well as crystal lattice 
distortions. 



Results and discussion 



If the JT contribution from lattice energy (Eq. ^ will be removed {\Ve\ = in Eq. |0J), the structure, obtained 
via minimization of energy, would have space group R3c (tabl. EJ. Adding JT term to the crystal energy we will get 
monoclinic phase with P2i/h symmetry. Structural phase transition from i?3c to P2i/b structure in rhombohedral 
LaA'InOs was observed We suppose that in the high temperature rhombohedral phase JT effect remains, but it is of 
cooperative dynamical type. We performed the calculation of energy dependence on the values of Eg vibrational mode 
components qg and q^. For each fixed values of these parameters the crystal energy was optimized with respect to the 
all lattice parameters which do not breaks new. Pi, symmetry. Because of fixed values of qg and qe parameters there 
was some constrains on the optimization parameters: a) the projection of ions coordinates onto qg, displacements 
were constant; b) the lattice constants were constrained to be like in i?3c lattice, i.e. the lattice was allowed only to 
change its volume and elongate along cubic [111] direction (such constrains do not allow for the uniform JT distortion) 

The calculated dependence could be described well with following relationship (energy per primitive cell) 

E=^ {ql + ql) - 2\V,\^ql + ql + A {Sq^e ql) (9) 

Where k = 12 eV/A^, \Ve\ = 1.3 eV/A and A = 0.005 eV/A^ Because of small value of A system will perform "slow 
rotation" at the bottom of adiabatical potential. The trajectory of motion will be approximately described by the 
equation \/ {qg + q^) = p = 2\Ve\/k. The value of p is much more grater then the oxygen ions vibrational motions 
amplitude. We suppose that the frequency of motion at the bottom of adiabatic potential is slower then frequencies 
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of optical phonons. We have to average over the trajectory of "motion" in the qs, qe space to get observable physical 
values. The RS process is "fast" process, but the total time of observation is much grater than process itself. That is 
why RS process "feels" each low symmetry structure that corresponds to different values of but the total picture 
is averaged on the whole trajectory of motion. In such case we have to average phonon frequencies and their RS 
intensities over the different $ values. The dynamic matrix is calculated and diagonalized for each value of $. As 
the result the averaged phonon frequencies can be got by the averaging over the $ dependence (table IIIip . The 
dynamic breaking of symmetry could give rise to the appearance of new RS lines. For verification of such assumption 
we calculated the dependence of intensity of highest vibration (the symmetry in i?3c structure) on the value of 
\Ve\- For each value of \Ve\ we calculated the p which corresponds to the energy minimum (the bottom of adiabatic 
potential). For this value of p the phonon frequencies and RS lines intensities were averaged over all values of <& 
(Fig. 13. 

Because of low symmetry at each point of trajectory of averaging the phonon frequencies of Eg symmetry are split 
(T ab, jmil . However, the value of gap is rather small to see distinct lines and to give rise for observed in experiment 
III lilt il4j broadening. The broadening of lines could have the vibronic origin, which is integral feature of systems 
with dynamical JT effect 0. 
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TABLE L The pair potential model parameters (all in atom units) 



A 85.623 32.689 44.862 X 9 8.265 4.245 

B 1.276 1.671 1.0592 Y -6 -6 -6 

C 250.0 43.56 21.729 k 40.0 50.76 3.99 

D 1.650 1.575 1.320 

AO 4.513 



TABLE II: The structure of LaMnOs in rhombohedral phase 





a 


c 


X 


y 


z 


Calc." 


5.809 


13.894 


0.455 
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Aver. Calc* 


5.837 


13.911 


0.451 





1/4 


Exp. Ref.^ 


5.529 


13.335 


0.447 
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"The crystal lattice parameters were calculated via energy optimization having | Ve | =0 

'The crystal lattice parameters were obtained via averaging of optimized structure having IVel = 1-29 eV/A 



TABLE III: RS lines frequencies 



Mode 



LaAlOs 



LaAInOa" 



LaMnOi" 



LaMnO-/- 



LaMnO/ 



Eg 

Alg 
A2g 
Ea 



12s 



33 

180 
123 
152 
487 



39 

144 
172 

223 
261 
312 

468 

580 



179 
236 
329 
515 
640 



217 

323 

427 
497 
618 



61 

74 

142 

168 

177 

235 

261 

314 

322 

e 

478 

584 



"The experimental values are taken from Ref. Q An assignment is carried out according to our intensities calculations. 
''The calculated frequencies without JT term in c rvstal energy and dynamical matrix {R3c space group) 
■^The experimental values are taken from Ref. Il4l 
■^An averaged on the AP bottom RS frequencies 

■^This mode corresponds to the motion along the bottom of adiabatical potential and to got frequency of it will require separate consid- 
eration 
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FIG. 2: The dependence of intensity of highest A2g (in rhombohedral phase notation) mode on the value of | Ve|. The frequency 
and intensity were averaged on "1> for the value of p, which corresponds to the energy minimum. For each value of <1> the ions 
were "relaxed" according to scheme described in text 



